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ABSTRACT 

The recent discovery of a neutron star accretor in the ultra-luminous X-ray source M82 X-2 challenges our 
understanding of high-mass X-ray binary formation and evolution. By combining binary population synthesis 
and detailed mass-transfer models, however, we show that the binary parameters of M82 X-2 are not surprising 
provided non-conservative mass transfer is allowed. Specifically, the donor-mass lower limit and orbital period 
measured for M82 X-2 lie near the most probable values predicted by population synthesis models, and systems 
such as M82 X-2 should exist in approximately 13% of the galaxies with a star-formation history similar to 
M82. We conclude that the binary system that formed M82 X-2 is most likely less than 50 Myr old and contains 
a donor star which had an initial mass of approximately 8-10 Mq, while the NS’s progenitor star had an initial 
mass in the 8-25 M© range. The donor star still currently resides on the main sequence, and is capable of 
continued MT on the thermal timescale, while in the ultra-luminous X-ray regime, for as long as 400,000 
years. 

Subject headings: stars: binaries — stars: neutron — X-rays: binaries 


1. INTRODUCTION 

Ultra-luminous X-ray Sources (ULX) are among the most 
extreme phases of binary evolution, characterized by X-ra y 
luminosities exceeding 10^^ erg s“^ (iGladstone et al.l[200^ . 
Since these systems exceed the Eddington luminosity for 
compact objects formed via stellar evolution, models must 
invoke one of two mechanisms in order to account for 
their high luminosities. Under the assumption of spherical, 
Eddington-limited accretion, one class of models employs 
intermediat e-mass black holes (BHs) w i th masses exceeding 
^ 100 Mr^ (IColbert & Mushotzkvil 199^: Ivan der Marelll2064 

iMadhusudhan et al.ll2008l: iFeng & SorialbOl ih . The second 

method employs some combination of thi n accretion disks 
(IShakura & Sunvae^l 197 3l: lBegelmanl[2002h . whose luminos¬ 
ity can exceed the Eddington limit without being disrupt ed by 
radiat ion pressure, and/or anisotropic X-Ray emission (iKingl 
1200^ . in order to reproduce the luminosity of observed sys¬ 
tems without exceeding the Eddington limit. The formation 
rate of ULXs has been accounted for in models examining the 
formation rate of th e most massive BHs from single stars (e.g. 
iMapelli et al.ll201Q|) and with the inco rporation of binary evo¬ 
lution effects te.g. lLinden et al.ll201Qh . 

Recently, NuSTAR discovered a pulsar spatially coinci¬ 
dent with the location of the ULX M82 X-2, indicating 
that the compact object in this system is a neutron star 
(NS; iBachetti et al.l l2014l) . Assuming a typical NS mass 
Mns ^ 1.4M0, the measured mass function and orbital pe¬ 
riod (2.52 days) indicate that the donor star has a mass M 2 > 
5.2 Mq and a radius R 2 ^ 7 R©. The X-ray luminosity of M82 
X-2 is observed to be Lx(0.5 - 10 keV) = 6.6 x 10^^ erg s“\ 
indicating that the system has an accretion rate that exceeds 
the Eddington limit by at least a factor of 30. 

The identification of M82 X-2 as a NS UL X has a lready 
generated a variety of new theoretical insights. iTongl (120141) 
finds that the spin-up behavior of the pulsar is consistent 


with a low- magnetic field (^ 1TG ) magnetar, while lEksi et al.l 
(120141) and iDairOsso et al.l (120141) use the physical properties 
of the system to imply la rger magnetic fields ^10-100 TG. 
iKluzniak & Lasotal (1201 5h note that the high MT rates would 
spin up the NS to millisecond periods within 10^ yr, open¬ 
ing the possibility that observable high-mass X-ra y binaries 
may b e found with a millisecond pulsar accretor. iLyutikovI 
(120141) suggests that the extremely super-Eddington accretion 
rates require a new MT regime, where an optically thick ac¬ 
cretion curtain shields the interior gas from the outgoing X- 
ray fiux. Furthermore, this analysis suggests that the equi¬ 
librium, where the Alfven radius matches the co-rotation ra- 
dius, indicates a inagnetic field of 14 TG. An analysis by 
IChristodoulou et all (120141) disagrees, arguing that the lumi¬ 
nosity of M82 X-2 can be explained just with mild geometri¬ 
cal beaming. 

Apart from the specific emission mechanism and the mag¬ 
nitude and influence of the NS’s magnetic field, the measured 
orbital period and the lower limits on the donor mass and ra¬ 
dius are both telling and puzzling at first glance: (1) the high 
inferred accretion rate onto the NS requires that the donor is in 
Roche-lobe overflow (RLO), as wind-fed X-ray binaries with 
NS accretors are expected to have X -ray luminosities or ders 
of magnitude below the ULX range (iLinden et al.ll2010l) . (2) 
the donor star must be hydrogen rich, as a helium star with 
mass > 5 Mq cannot fill its Roche lobe in a 2.5 day period or¬ 
bit, independent of its evolutionary stage, (3) the evolutionary 
mechanism must have either a long lifetime, or a high for¬ 
mation rate, in order for this one system to exist in the local 
universe. 

Still, the picture of a donor more massive than ~ 5M0 
transferring mass onto a NS is particularly puzzling. While 
RLO X-ray binaries with NS accretors are a common output 
of population syn thesis models, they typically invo lve lower- 
mass donors ('e.g. lFragos et al.ll200Ml2009aLl2013allE) . This is 
due, primarily, to the stability criteria for RLO from a mas- 
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sive donor star onto a lower-mass accretor. Early polytropic 
and semi-analytic models, with the assumption of conserva¬ 
tive MT, demonstrated that systems with a donor more than 
times as massive as the compact object would quickly 
enter a delayed dynamical ins tability (iHjellming & Webbinld 
119871: llvanova & TaamI 120041) . producing a stellar merger. 
However, more recent studies of the stability of binary MT 
phases, using detailed binary evolution codes, showed that 
some of the approximations made in earlier studies, such as 
fully adiabatic mass-loss and the strict enforcement of hy¬ 
drostatic equilibrium, underestimate the maximum mass-ratio 
of a binary in which dynamically stable MT can occur (e.g. 
llvanova et al.1120131: iPavlovskii & Ivanovall20f^ . A crucial 
factor in determining this maximum mass ratio for stabil¬ 
ity of non-conservative MT phases is the specific angular 
momentum of the material ejected from the binary. Sev¬ 
eral different models have been develo ped for different physi¬ 
cally motivated mass-loss models (e.g. Ivan den Heuvellll994 
iKing & Ritterlll99^lTauris & van den Heuvelll2006l). 

In this Letter, we investigate the formation of ULX systems 
with NS accretors and derive constraints on the evolutionary 
history and current properties of M82 X-2. We first utilize the 
Binary Stellar Evolution (BSE) code (iHurlev et al.l 120021) in 
order to study the initial binary properties that can produce a 
MT event with characteristics similar to M82 X-2. We termi¬ 
nate simulations with this code at the onset of RLO and per¬ 
form the MT calculations with the de tailed stellar evolution 
code MESA dPaxton et al .11201 lLl2?m . in order to accurately 
predict the duration and accretion rate of the MT phase. 

2. THE HOST GALAXY: M82 

M82 is one of the nearest dusty starburst galaxies. It 
has sustai ned vigorous star f ormation in the past 50 to 
100 Myr (iBarker et al.l l2008l) . probably triggered by its 
anterior interaction with M81. Based on publicly avail¬ 
able ph otometry from IRA S, Spitzer-MIFS, and Herschel- 
SPIRE dRoiissel et alJl2(mih . we estimate the IR-derived star- 
formation rate (SER) to be of at least 5.5Mpvr~^ d uring 
this period, in agreement with iKennicutt et air(l2003l) . We 
should note, however, that preliminary spectral energy dis- 
tri bution fitting results , following the methodology described 
in ISklias et HT (120141) . favor a declining star-formation ac¬ 
tivity in M82 over a constant SER. Along the same lines, 
earlier, spatially resolved, studies infer SERs as high as 
^ 30Mr7^vr~\ depending on various IME assumptions (e.g. 
lEorster Schreiber et all 120031) . In the rest of this paper, we 
adopt a constant SER of 5 M^;^ yr"^, over the last 100 M yr, and 
a metallicity of Z = 0.02 (lEorster Schreiber et al.ll2001l) for all 
our simulations. 

3. NUMERICAL MODELS 
3.1. Binary Population Synthesis with BSE 

In order to determine the most probable initial parameters 
for M82 X-2, we calculated the evolution of 10^ isolatecQ 
massive binaries, restricting our modeling to systems with 
a primary mass large enough t o form a NS. We utilize the 
BSE population synthesis code (iHurlev et al.ll2002l) . modified 
to include the su ite of stellar wind presc riptions for massive 
stars described in lBelczynski et al.l (I201QI) . and the fitting for¬ 
mulae f or the binding energy of the envelopes of stars de¬ 
rived bv iLoveridge et al.l (120111) . We followed binary systems 

^ Eormation channels through dynamical interactions in dense stellar sys¬ 
tems are not taken into account 


through any binary interaction prior to the primary’s super¬ 
nova, through the supernova explosion which formed the NS, 
and the post-supernova detached evolution. We halted the 
code at the moment when the secondary star fills its Roche 
lobe, and recorded the state of the binary to use as an input to 
our stellar evolution and MT code^ _ 

In our modeling, we assume a lKroupal (120011) initial mass- 
function, a flat distribution of initial binary mass-ratios, a 
logarithmically fiat distribution of initial orbital separations, 
a Maxwellian distrib ution of NS supern ova kick magnitudes 
with a = 265kms“^ (iHobbs et al.l 120051) . and a common en¬ 
velope efficiency of acE = 0.5. These parameter values have 
been shown to produce synthetic populations of X-ray bina¬ 
ries that are in good agreement with observations of X-ray lu¬ 
minosity functions an d X-ray scaling relations of extragalac- 
tic populations, (e.g. lEragos et al.l 120081: iLinden et al.l [20091: 

lEragos et'n]l2013al: iTzanavaris et al.ll2013l) . We repeated all 

calculations in this paper utilizing common envelope efficien¬ 
cies of acE = 0.2 and acE = 1-0, finding only a small change 
in the resulting distributions for the orbital periods and donor 
masses of possible M82 X-2 progenitors. The qualitative in¬ 
sensitivity of our models to the common envelope efficiency 
is explained by the high mass of the secondary star (> 5 M 0 ) 
which guaranties that the orbital energy of the binary is usu¬ 
ally greater than the binding energy of the NS’s progenitor 
envelope. 

3.2. Detailed Mass-Transfer Calculations with MESA 

We use version 7184 of the stellar evolution code MESA 
(iPaxton et al.l 1201 iL 1201 3l) in order to calculate grids of ^ 
2,700 MT sequences for NS X-ray binaries undergoing MT. 
We assume an initial NS mass of 1 . 4 M 0 , and cover the 
parameter space of possible orbital periods and companion 
masses at the onset of RLO. Specifically we consider initial 
donor masses (M 2 ) between I.OM 0 and 5 O.OM 0 and initial 
binary orbital periods between 1.0 and 300.0 da y. A higher 
initial NS mass (e.g. 2.0M 0 ), as suggested bv lTauris et al.l 
(I 2 OIII) for the formation of the binary millisecond pulsar PSR 
J1614-2230, will yield in our case qualitatively identical re¬ 
sults, and is not considered in this work as it would increase 
significantly the computational cost of our simulations. 

Eor our MT calculations we employ the implicit MT 
scheme and the hydrodynamic solver of MESA, in order to 
properly resolve the onset of a dynamically unstable MT 
phase. We terminate the MT sequences when any of the fol¬ 
lowing criteria are met: (i) the mass of the donor star drops 
below 1 M 0 , (ii) the age of the system exceeds the age of the 
Universe (13.7 Gyr), (in) the whole envelope of the donor star 
is removed, or (iv) the timestep of the hydrodynamic solver 
drops below 10“^ s, while the MT rate exceeds lO“^M 0 yr“\ 
and no solution can be found such that the radius of the donor 
is less or equal to the radius of its Roche lobe. We consider 
the latter as the onset of a common envelope. 

The secular evolution of the orbit of a mass-transferring bi¬ 
nary depends on the accretion efficiency, and on the specific 
angular momentum tha t is los t from the system. Here, fol¬ 
lowing [yanden^eu^ (II 994 I) . we define a parameter a de¬ 
noting the fraction of the mass lost from the donor (M 2 ) that 
is ejected from the vicinity of the donor star in the form of 
“fast wind”, i.e. with specific angular momentum equal to 
that of the donor star. The remainder of the mass ((l-a)M 2 ) 
is funneled through the first Lagrangian point to the accretion 
disk around the NS. We limit the accretion rate onto the NS 
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Figure 1. Left panel: The two-dimensional distribution of binary orbital periods versus donor masses of binaries that contain a NS and a secondary non¬ 
degenerate star that just filled its Roche lobe, calculated from 10^ isolated binaries evolved with BSE. Right panel: The gray shaded area shows the parameter 
space in the orbital period-donor star mass plane that is covered by our grid of detailed MT sequences calculated with MESA, assuming a MT efficiency parameter 
OL = 0.9. 

ceeds 


to 100 times the Eddington limit We note here that 

the exact value of this upper limit, within a factor of a few, 
does not play a significant role in the evolution of a mass- 
transfer sequence. The fraction of (1 -ct)M 2 that is in excess 
of lOOM^jj is again ejected from the system in the form of 
“fast winds”. However, we assume that this mass is ejected 
from the vicinity of the NS and has the specific angular mo¬ 
mentum of the accretor. Since these X-ray binaries have a 
donor star that is more massive than the NS, mass lost from 
the vicinity of the donor star has a smaller angular momentum 
than mass lost from the vicinity of the accretor. Hence, a value 
(T > 0 would results in a reduced rate of angular momentum 
loss from the orbit, and would therefore increase the stability 
of the MT phase. We ran our grid of MT sequences for 4 dif¬ 
ferent values of the parameter a(a = 0.0,0.5,0.9, 0.99) in or¬ 
der to study the effect that it has on the formation of potential 
NS ULXs. The right panel of Figure [T] shows the parameter 
space in the orbital period-donor star’s mass that is covered 
by our grid of detailed MT sequences, assuming an initial NS 
mass of 1.4M0 and a MT efficiency parameter ol = 0.9 

4. RESULTS 

In Figure \T\ (left) we show the distribution of orbital pe¬ 
riods and donor masses, at the beginning of the MT phase, 
as produced by the BSE code. Intriguingly, we note that the 
peak of this distribution is centered around systems with a 
donor mass of 5-10 Mq and an orbital period between 1.6 
-3.2 days at the onset of RLO. This distribution of binaries 
is thus highly compatible with the observation of M82 X-2. 
The population appears to change in character at ^13 Mq; 
systems to the left of this boundary experience a CE phase 
before NS formation, while systems to the left evolve through 
stable MT between the two initial massive stars. Similarly, 
there is a “valley” around orbital periods of 3-10 days, which 
separates systems with main-sequence donors from systems 
with evolved donors. However, in order to estimate the prob¬ 
ability that any of these systems is currently a ULX, we need 
to convolve this distribution function of binaries at the onset 
of RLO with the detailed MESA calculations, so that we de¬ 
termine the duration of ULX activity in each binary. Eor our 
analysis using the MESA code, we define as ULX any X-ray 
binary in which the NS is accreting material at a rate that ex- 

^ For the calculation of MEdd a radiative efficiency of 10% is assumed. 


In Eigures [2] and [3] we show the resulting probability dis¬ 
tribution of having a ULX with a NS accretor in the parame¬ 
ter space of the observed characteristics of M82 X-2 (donor 
mass and orbital period). We show results for four differ¬ 
ent age ranges for the ULX, 0-25 Myr, 25 -50 Myr, 50 - 
75 Myr, and 75 - 100 Myr, and four different values of non¬ 
conservative MT: a = 0.99, 0.9, 0.5 and 0.0. The color bar 
depicts the relative probability of a NS ULX existing today 
with the given orbital parameters, and the expected number 
of systems in M82 is listed as < A^ns-ulx >• This probability 
is calculated by assigning to each one of our MT sequences 
a weight proportional to the formation rate of NS binaries 
reaching RLO at periods and donor masses close to that of 
the MT sequence, as calculated by our BSE population syn¬ 
thesis models. In addition, we weight each MT sequence by 
the time that the sequence spends at each period and donor 
mass, while the MT rate exceeds lOMEdd- 

In each panel of Eigures [2] and [3l we show white lines de¬ 
noting the observed orbital period and minimum donor mass 
of M82 X-2. We note two key observations: (1) the observed 
parameters of M82 X-2 are well matched to the most probable 
parameters of a ULX with a NS accretor, (2) the probability 
of such a system in M82 existing is relatively high, between 
~ 0.03 and ~ 0.13 (i.e. we expect to have a NS ULX in every 
one out of approximately 8 galaxies, with properties similar 
to M82). Eurthermore, comparing the observed properties of 
M82 X-2 to the different panels of Eigures [2] and [3l we can 
put limits on both the current age of the system and the accre¬ 
tion efficiency of the MT phase. Accretion efficiencies below 
1% (a > 0.99) can be excluded, as the accretion rate becomes 
too low to produce a ULX (see the very low expected number 
of ULX in the all panels of the first row of Eigure[2]). Simi¬ 
larly, accretion efficiencies above 50% (a < 0.5) predict that 
ULXs would be formed at significantly shorter periods than 
observed in M82 X-2. The age of the population plays also 
role, as older systems tend to have shorted periods and lower 
mass companions. Based on this comparison, a system like 
M82 X-2, would likely have an age < 75 Myr, and an accre¬ 
tion efficiency below 10%. 

We can derive more robust constraints on both the current 
properties of M82 X-2 and the properti es of its progenitor bi¬ 
nar y by following an analysis sim ilar to lEragos et al.l (l2009bh 
and lEragos & McClintocla (120151) . Eor each MT sequence in 
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Figure 2. Relative likelihood of a NS ULX existing today at a given orbital period and donor mass. Each column of the hgure corresponds to a different age 
range of the ULXs, while each row corresponds to a different choice of the accretion parameter a (top: a = 0.99, bottom: a = 0.9). 
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Figure 3. Same as Figure|2] but for a. = 0.5 (top) and a. = 0.0 (bottom). 

our grid, we examine all systems that simultaneously fit the 
observational constraints on M82 X-2, specifically an orbital 
period of 2.52 day, a donor mass above 5.2 Mq and an accre¬ 
tion rate onto the NS higher than lOM^jj. In this part of the 
analysis we only consider one of our MT grids, with ct = 0.9, 
which predicts the highest expected number of NS ULXs and 
the closest-matching properties for M82 X-2. In Figure [4] we 
show the probability density functions of current binary pa¬ 
rameters (donor mass and age, and ULX duration) for sys¬ 
tems with the M82 X-2 measured characteristics. We also use 
the results from our parameter-space exploration in BSE in 
order to derive self-consistent constraints on the initial binary 
parameters of M82 X-2’s progenitor (component masses and 
orbital period). 


5. DISCUSSION 

By combining a parameter space exploration using the BSE 
population synthesis code with a detailed treatment of MT us¬ 
ing the MESA code, we have shown that the properties of 
M82 X-2 are well explained by current models of binary pop¬ 
ulation synthesis. Specifically, our models reveal the follow¬ 
ing details about the formation of M82 X-2: 

• The donor star is hydrogen rich, as a helium star of mass 
> 5.2M0 cannot h\\ the Roche lobe of a ^ 2.5 day or¬ 
bit, even at the late stages of the carbon/oxygen burning 
phase. Eurthermore, the donor star is in RLO, as wind- 
fed NS binaries are expected to have X-ray luminos ities 
below the Eddington limit fe. g. iLinden et al .Il20 1 Ql) . 
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Figure 4. Probability density functions of the current mass and age of the companion star in M82 X-2, the initial component masses and orbital period of M82 
X-2’s progenitor binary, the orbital period of the binary at the onset of RLO, as well as the duration of the ULX phase and the current mass-loss rate of the 
companion star. We take into account observational selection effects by weighting all probability density functions with the total time that each modeled system 
has properties similar to M82 X-2, i.e. M 2 > 5.2Mq, Maccvetod > lO^Edd. and 2 < < 5day. Finally, in the bottom right panel we show for comparison the 

MT rate corresponding to the Kelvin-Helmholtz thermal timescale of an 8 M 0 star that is just hlling its Roche lobe in a 2.5 day orbit around a 8 M 0 NS. 


• The orbital period of any ULX with a NS accretor is 
most likely to be observed with an orbital period be¬ 
tween 1-3 days, and a donor mass between 3-8 Mq. 
This places the observed properties of M82 X-2 near 
the peak of the likelihood distribution for this class of 
systems. 

• The MT is highly non-conservative and happens on 
the thermal timescale. Conservative MT leads to 
dynamical instability (p ossibly delayed) as discussed 
in the literature (e.g . iHjellming & Webbing 119871: 
llvanova & TaamI120041) . The accretion efficiency must 
be < 0.1, consistent with local ULX analogues (e.g. 
iNeilsen & Leell200^ iPonti et al.ll2012h . and mass lost 
should be ejected from the system in the form of 
“fast winds” from the vicinity of the donor, i.e. it 
should carry relatively little angular momentum in or¬ 
der to provide stability. A good example of this is 
the short-period high- mass X-ray binary Cygnus X-3 
(Ivan den Heuvell[l994l) . 

• Assuming an accretion efficiency of 10%, we esti¬ 
mate that the number of NS ULXs per unit of SFR is 
Nns-ulx/SFR = 0.027 MQyr. This number is an order 
of magnitude lower compared to predictions in the lit¬ 
erature for the formation rate of ULXs with BH accre- 
tors, which sugges t a Nbh- ULX /SFR ^ Q;2~ L OM^Ijt 
at sol ar metallicity (iRappaport et al .1120051: iLinden et al.l 
I 2 OIOI) . We stress here that for a NS ULX to be identi¬ 
fied as such, the requirement of a highly magnetised 
NS, that allows the production of X-ray pulses, should 
be taken into account. The modelling of the NS’s mag¬ 
netic field is outside the scope of this work. 
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